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An explanation for the soft X-ray excess in AGN 
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ABSTRACT 

We present a large sample of type 1 active galactic nuclei (AGN) spectra taken with 
XMM-Newton, and fit them with both the conventional model (a power law and black 
body) and the relativistically-blurred photoionized disc reflection model of Ross & 
Fabian (2005). We find the disc reflection model is a better fit. The disc reflection 
model successfully reproduces the continuum shape, including the soft excess, of all 
the sources. The model also reproduces many features that would conventionally be 
interpreted as absorption edges. We are able to use the model to infer the properties 
of the sources, specifically that the majority of black holes in the sample are strongly 
rotating, and that there is a deficit in sources with an inclination > 70°. We conclude 
that the disc reflection model is an important tool in the study of AGN X-ray spectra. 

Key words: accretion, accretion discs - galaxies: active - X-rays: galaxies. 



1 INTRODUCTION 



Active galactic nuclei are classified into two types: type 
1 AGN have both broad and narrow components of their 
optical emission lines, type 2 AGN show only narrow lines. 
The narrow lines are emitted from the outer regions of the 
source, the broad lines from matter in the inner regions 
where Doppler broadening increases the width of the lines. 
The unified model of AGN explains the two types as an 
inclination effect; Seyfert 2 galaxies are observed at high 
inclinations and have the central regions obscured by a 
torus, Seyfert 1 galaxies are seen at low inclinations where 
the line of sight passes over the torus. Therefore, we expect 
that our sample of type f AGN will allow us to investigate 
the central regions of these systems. 

The soft excess is a major component of the spectra of 
many AGN in the X-ray band and is present in every source 
in this survey, see Figure First noted in Arnaud et al. 
(1985), it is extra emission below ~1.5 keV than expected 
from extrapolating the power law spectrum observed at 
higher energies. The soft excess is usually well fit with a 
black body which has a roughly constant temperature of 0.1 
- 0.2 keV over several decades in AGN mass (Walter & Fink 
1993, Czerny et al. 2003, Gierlifiski & Done 2004). If it is 
thermal, this temperature is much too high to be explained 
by the standard accretion disc model of Shakura & Sunyaev 
(1973), although it could be explained by a slim accretion 
disc in which the temperature is raised by photon trapping 
(e.g. Abramowicz et al. 1998; Mineshige et al. 2000) in 



E-mail: jc@ast.cam.ac.uk 



which case the accretion is super-Eddington (Tanaka, Boiler 
& Gallo 2005), or by Comptonisation of EUV accretion 
disc photons (e.g. Porquet et al. 2004). Atomic physics 
could trivially explain the constant temperature of the 
component, e.g. strong, relativistically-blurred absorption 
from a disc wind (Gierlihski & Done 2004) or photoionized 
emission blurred relativistically by motion in an accretion 
disc (Ballantyne, Iwasawa & Fabian 2001; this paper). 
Relativistically-blurred photoionized iron emission has been 
observed in several AGN: a clear relativistic profile in a few 
sources (e.g. MCG -6-30-15, Tanaka et al. 1995; Vaughan & 
Fabian 2003), and extended red wings to iron lines in many 
sources (e.g. Nandra et al. 1997). This concept is extended 
in disc reflection models which reproduce the emission 
expected from a photoionized accretion disc around a black 
hole (e.g. Ross & Fabian 1993; Ballantyne, Ross & Fabian 
2001), models which have been used with some success in 
e.g. Ballantyne, Iwasawa & Fabian (2001), Fabian et al. 
(2002), Ballantyne, Vaughan & Fabian (2003). The work 
presented here is the largest study to date, and is focused 
on fitting many objects with a uniform analysis rather 
than a detailed study of one or a few sources. We uses 
the latest model from Ross & Fabian (2005), in which a 
semi-infinite slab of optically thick cold gas of constant 
density is illuminated by a power law, producing a Compton 
component and fluorescence lines from the ionized species 
in the gas, the illuminating and reflected components are 
added together (with the amount of each free to vary) and 
the total emission convolved with a Laor (1991) line profile 
to simulate the blurring from a relativistic accretion disc. 
The Laor profile assumes a flat accretion disc with a power 
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law emissivity and sharp edges at the inner and outer radius 
around a maximally-rotating black hole. The Laor profile 
has as free parameters the emissivity index, the inner and 
outer radius (in gravitational radii; the blurring process 
does not depend on black hole mass) and the inclination 
of the disc to the line of sight. The reflection spectra 
allow fitting to the iron abundance of the gas, the spectral 
index of the illuminating power law, and the ionization 
parameter (a measure of the ratio of the energy density in 
the illuminating radiation to the atomic number density 
in the illuminated gas). Examples of the model are given 
in Figs. |21 and 1131 This model is physically motived, and 
therefore may prove useful in understanding these sources. 
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2 THE SOURCES 

In this paper we use publicly available archival XMM- 
Newton data on 22 type 1 AGN from the Palomar-Green 
(PG) sample and a selection of 12 other Seyfert 1 galaxies 
with high-quality observations available. Where multiple 
observations are publicly available we use the longest where 
the pn took data. A list of the sources with their properties 
and the observation IDs used is given in Table Most of 
these sources are the subject of many papers and have been 
observed with several different X-ray instruments, however, 
we confine ourselves to XMM-Newton data and perform a 
standard analysis across all of them. 



3 DATA REDUCTION 

We obtained the Observation Data Files (ODFs) from the 
XMM-Newton public archive and reduced them in the stan- 
dard way using SAS 6.0 to produce event lists. We extracted 
spectra and light-curves using circular source-centred re- 
gions 40 arcsec in size for the pn and 60 arcsec in size for 
the MOS, using a smaller extraction region where appropri- 
ate due to chip gaps, etc. We used the standard valid event 
patterns of - 4 (singles and doubles) for pn and 0-12 
(singles to quadruples) for the MOS. Background spectra and 
light-curves were similarly created using regions away from 
any sources, possible out-of-time event trails, and chip gaps. 
Where the count rate of background flares was greater than 
5 per cent of the source count rate a Good Time Interval 
(GTI) file was used to exclude those events, and even more 
conservative background filtering was adopted where possi- 
ble. We used the SAS tasks rmfgen and arfgen to create 
the response matrices. 

As the MOS cameras are not accurately calibrated against the 
pn (Kirsch et al. 2004) only one of the two systems can be 
used without introducing complications to the fits. We chose 
to use pn data as it has a larger effective area and spectral 
range. We did not reduce RGS or 0M data. 
The SAS task epatplot was used on the extracted spec- 
tra to test for the presence of pile up. All but five of the 
observations were found to be pile-up free. The piled up ob- 
servations were re-extracted using an annular region to ex- 
clude the piled-up centre of the observation. Annulus sizes 
used were 11, 8, 8, 16 and 5 arcsec for PG 0003+199, PG 
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Figure 2. The rclativistically-blurred photoionized disc reflec- 
tion model fit and residuals to PG 1501+106 (top) and a plot of 
the model used showing the illuminating power law and reflection 
components (bottom). The dashed line is the reflection compo- 
nent, the dash-dotted line is the power law, and the solid line is 
the total model. The model is shown over an extended energy 
range. The curvature at low energies is due to cold absorption 
local to our Galaxy. See Table 131 for fit parameters. 



0844+349, PG 1244+026, PG1440+356 and RE J1034+396, 
respectively. 

We performed the spectral analysis using xspec vll.3 (Ar- 
naud 1996), taking the limits in accurate calibration of the 
pn data as 0.3 - 12.0 keV. We grouped the spectra so each 
bin would include 20 or more source counts so that \ 2 statis- 
tics would be applicable. All quoted errors are 90 per cent 
limits on one parameter (A^ 2 = 2.706). 
We produced background-subtracted light-curves and hard- 
ness ratios for the sources. We found that 12 of the 34 sources 
display some variability in their hardness ratios over the 
course of the observations. In the interest of using the high- 
est quality data available, we analysed only spectra inte- 
grated over the entire observation. This strategy may have 
introduced some error, but a detailed analysis of spectral 
variability of a large number of sources is beyond the scope 
of this paper. 



4 ANALYSIS AND RESULTS 

We first fit the spectra with a power law and redshifted (lo- 
cal to the source) black body, modified by absorption from 
the cold gas in our Galaxy, with the amount of cold absorp- 
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Figure 1. The soft excesses of all the sources in our sample. The sources were fit as in Table 1^1 then the black body component was 
removed and the residuals were plotted. Absorption edges were left in the model to show the form of the soft excess more clearly. The 
top three sources, in descending order of flux at 0.5 keV, are ARK 564, NGC 4051 and TON S180; note the unusually extended shape 
of the TON S180 soft excess. 



tion fixed at the value given in Table □ We then allow for 
variable amounts of redshifted cold absorption (i.e. due to 
cold gas in the source galaxy) and warm absorption from 
ionized gas in the line of sight, and finally narrow iron emis- 
sion due to e.g. reflection from a cold torus 1 . We then fit 
the same data with the relativistically-blurred photoionized 
disc reflection model 2 , which we likewise check for the pres- 
ence of redshifted cold and warm absorption, and narrow 
iron emission. We also investigate the possibility that the 
illuminating power law component is absent. Where adding 
a component produces an improvement in x 2 of < 2.7 per 
lost degree of freedom we do not include that component in 
our final fit. 

Specifically, to check for warm absorption we use two red- 
shifted absorption edges whose energies are initially fixed 
at 0.74 and 0.87 keV to model the O VII and O VIII K 



1 In XSPEC terminology our simple model is 
phabs*zphabs*zedge*zedge* (power law+zbbody+zgauss) 

2 Expressed in xspec commands our relativistically- 
blurred ionized disc reflection model is 
phabs*zphabs*zedge*zedge* (kdblur (powerlaw+ref lion)+zgauss) , 
where kdblur is the relativistic convolution and reflion is a 
table model of the ionized reflection. Note that reflion includes 
redshift as a parameter. 



absorption edges, as oxygen is the most abundant element 
in the soft X-ray band. We also allow the edge energy to 
vary to model in- or outflows, rejecting any fits outside 
0.45 - 1.1 keV as being unlikely to originate from O VII 
or O VIII. Further to this, given the recent controversy over 
whether AGN absorption edges are local to the sources or 
due to absorption from gas near our Galaxy (e.g. McKer- 
nan, Yaqoob & Reynolds 2004) we also test the possibility 
of Galactic oxygen absorption with non-redshifted edges at 
0.74 and 0.87 keV (we do not allow these energies to vary 
as relativistic outflows from our Galaxy are unlikely). We 
check for narrow iron emission using a redshifted Gaussian 
(A(E) oc e -(E-E f/2a 2 ^ with line width a fixed at 5 eV ( go 

the width is determined by the instrument resolution since 
the instrumental width is 3> the intrinsic line width) and 
energy Eq initially fixed at 6.4 keV (the energy at which 
neutral iron fluoresces most strongly), then allowing the en- 
ergy to vary to account for the possibility of in- or outflows 
and ionized iron fluorescence. Since we are fitting the entire 
(0.3 - 12.0 keV) spectrum whilst we test for the narrow iron 
line we only report lines for which xspec calculates a non- 
zero minimum equivalent width, this avoids the possibility 
that a detected line is due to a curvature in the continuum 
which our model has not addressed. Our technique is a con- 
servative one and is likely to miss marginal lines, we only 
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Table 1. List of sources investigated. The Galactic absorption column is taken from the NH ftool (Dickey & Lockman 1990). 
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a Mass taken from Gierlifiski & Done (2004), who derive them from Boroson (2002). b Mass from Woo & Urry (2002). c Mass from 

Wang & Lu (2001). 



include lines which affect the overall goodness of fit. The 
comparison in goodness of fit between the two models will 
not be affected by any features that both models do not re- 
produce. A survey designed to detect fainter features (using 
the F-test), and including many of the same sources is avail- 
able in Porquet et al. (2004). 

We present the best black body and relativistically-blurred 
photoionized disc reflection model fits in Tables |5] and |3] re- 
spectively. A plot of our measured soft excess temperature is 
shown in Fig. and several example comparisons between 
the two different model fits are plotted in Figs. 0] and |H] 
Example vF v plots are shown in Fig- [HI 
We find cold absorption in only one source, PG 0050+124, 
and clear narrow iron lines in two, NGC 4051 and MRK 
766. The same sources show extra cold absorption local to 
the source and narrow iron lines with either model, with the 
exception of PG 1211+143 which only appears to have an 
iron line when the simple model is used. PG 1211+143 shows 
spectral complexity in the 4-8 keV region, see Section |K| 
and Fig. |S] We find redshifted warm absorption edges in 18 
sources using the simple model, compared to 7 sources us- 
ing the disc reflection model. We find non-redshifted warm 



absorption in only one source using the simple model, and 
three sources using the disc reflection model. We do not de- 
tect a power law component in 11 of the 34 sources when 
using the disc reflection model. 

Note that the disc reflection model has some limits built in. 
The Laor profile has an inner disc radius that cannot be < 
1.235 gravitational radii, for physical reasons discussed in 
Section|S] The profile also has a disc emissivity index which 
cannot be > 10, this is quite reasonable as an emissivity 
profile going as the tenth power of radius (n.b. where r is 
radius and A is area r~ 10 dA oc r~ 9 dr since dA/dr — 2%r 
for a disc) emits 90 per cent of its radiation between 1.235 
(the inner edge) and 1.65 gravitational radii (this may not 
equal the fraction that reaches the observer due to various 
relativistic effects, see Section The reflection model is a 
grid model which is restricted in parameter space to increase 
the speed with which it can be computed. The limits in our 
version of the model which the sources encounter in fitting 
are the constraint that iron abundance (Fe) must be < 10 
times solar and the ionization parameter (£) is restricted to 
1 < ^ < 10000 erg cm s _1 . 

It is also possible to create a disc reflection model based on 
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a non-rotating black hole, by using the diskline (Fabian et 
al. 1989) profile instead of the Laor profile for the blurring. 
This model is less successful at fitting the data; whilst it 
matches the general shape of the spectra it is less smooth 
than the Laor-blurred model and leaves residuals. We also 
fit using this diskline model, for all the sources the fit is 
worse than the Laor based model and is completely inad- 
equate (xt > 2) in several cases. We report sources where 
the fit is not significantly worse in Table |S] 



4.1 Reflection fraction and flux fraction 

Conventionally, the reflection fraction is the ratio of re- 
flected emission to power law emission. It is based on the 
assumption of an isotropic source of power law emission, 
a certain fraction of which is reflected from gas. The 
reflection fraction is then trivially related to the solid angle 
at the source of power-law emission which is covered by 
reflecting gas. An infinite plane of gas illuminated by an 
isotropic source would have a reflection fraction of 0.5. A 
source with no reflection has a reflection fraction of zero 
and a source with no power law emission has an infinite 
reflection fraction. Reflection fractions above one show that 
the AGN geometry is more complicated than a source and 
some clouds of gas, implying e.g. anisotropic emission or 
obscuration of the power law source. 

In the vicinity of a black hole, many of the fundamental 
assumptions used to calculate a reflection fraction and 
extract physical meaning from it are known to be violated. 
The matter in the disc is moving at relativistic speeds, 
so reflected radiation from the disc tends to be Doppler 
beamed in the plane of the disc, violating isotropy. In a 
non-relativistic situation all the energy that is incident on 
the disc is reprocessed and reemitted, but with Doppler 
boosting the reemitted radiation as detected by an observer 
may pick up more energy from the motions of the disc. 
Gravitational redshifting likewise makes comparing emitted 
and reflected energy non-trivial. Finally, space is curved, 
and light no longer propagates in straight lines - even if the 
source of power law emission is isotropic and not moving at 
relativistic speeds (i.e. with no beaming), the light bending 
effects will ensure that more radiation impacts on the disc 
than is detected by an observer, again violating isotropy 
(e.g. Miniutti & Fabian 2004). All of these effects may 
either increase or decrease the observed reflection fraction, 
and complicate the situation to such a degree that the 
conventional reflection fraction becomes physically mean- 
ingless. We therefore report only the purely observational 
quantity we call the flux fraction, which we define as the 
measured fraction of reflected flux over the total observed 
flux. 

We measure the flux fraction for each of the fits listed 
in Table [3] We use the 'flux' command from xspec to 
calculate the total flux of the model over the range 0.3 - 
12.0 keV, setting all warm and cold absorption and narrow 
Gaussian emission to zero so as to measure the flux emitted 
by the inner portion of the AGN, then set the power law 
component to zero and repeat the measurement to find 
the amount of reflected radiation. The flux fraction is the 
ratio of reflected emission to total emission, so a source 
with no observed reflection has a flux fraction of zero and 
a source with no observed power law emission has a flux 
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Figure 3. Soft excess black body temperature (kT) measured 
for the simple model plotted against black hole mass (M), see 
Section El and Table [TJ The soft excess is confined to a narrow 
temperature band over several decades in mass, as earlier noted 
by Gierlinski & Done (2004) and others. 



fraction of one. The previously mentioned infinite plane of 
gas illuminated by an isotropic source, with no relativistic 
effects, would have a flux fraction of approximately 0.5 
(not exactly 0.5 as reflection can shift photons outside 
the observed band). The flux fraction is dependent on the 
energy range it is measured over, so we use the same limits 
as our spectral fitting. The flux fractions are shown in 
Figure HJ 

We find that no source in our sample has a flux fraction be- 
low 0.25, that the flux fraction is roughly evenly distributed 
between 0.25 and 0.8, and that ~ one third of our sources 
have no power law component. This extreme reflection may 
be due to several effects, discussed in Section [6.71 We only 
measure one flux fraction above 0.8 and below 1.0 (with 11 
sources having a flux fraction of 1.0), this may be partly 
due to our fitting procedure; adding a power law component 
increases the number of degrees of freedom by one, so it 
must improve x 2 by > 2.7 for us to judge it significant and 
include it. Small power law components therefore tend not 
to be included, e.g. in our fit to PHL 1092 (a low quality 
observation where our best fit does not include a power law 
component) adding a power law component such that the 
flux fraction is 0.98 worsens the fit by a A\ 2 of 2, so a small 
power law component does not affect the fit enough to be 
significant. However, adding a power law component such 
that the flux fraction is 0.9 makes the fit worse by Ax 2 
of 11, enough to be clearly ruled out as a worse fit. The 
deficit in measured flux fractions between 0.8 and 0.9 is 
therefore not due to our fitting methods, but it may still be 
an artefact of our small sample size. It is also possible that 
there are two populations of AGN: those with flux fractions 
of 0.2 - 0.8 and those totally dominated by reflection. 
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Table 2. Simple fits to the sources. T is power law spectral index, kT is black body temperature, Edge Energy and r are the energy 
and optical depth of the two absorption edges. Energies given without associated errors are fixed at 0.74 or 0.87 keV, the energy of the 
O VII and O VIII K absorption edges. Blank spaces indicate no feature is present. 
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a The fit includes redshifted cold absorption with an equivalent hydrogen column of 3.8l°j' 2 x 10 2 ° cm 2 . b These fits include a narrow 
Gaussian line at 6.4 keV, equivalent width 80+.g° eV for PG1211+143, 109+.20 eV for NGC 4051 and 70 i30 eV for MRK ° 766 - 



5 NOTES ON INDIVIDUAL SOURCES 
5.1 PG 0050+124 

PG 0050+124 (also known as I Zw 001) has a very weak 
soft excess. When the simple fit is used without a black 
body component the fit is only slightly worse: xl °f 1-275 
for 930 degrees of freedom, compared to a xi of 1-270 for 
928 degrees of freedom. There is a Ax 2 of only 7 for 2 
fewer degrees of freedom when the black body component 
is added to the fit, so it is barely significant. 



5.2 TON S180 

TON S180 has an unusual soft excess which Vaughan et 
al. (2002) successfully modelled as a power law rather than 
a black body (see Figure 0. Neither of our models give 
an acceptable fit to the data, due to strong residuals in 



the pn spectrum below 0.5 keV (also present in the M0S 
spectrum, although of a slightly different shape). Vaughan 
et al. do not fit the pn spectrum below 0.5 keV, so their 
paper does not investigate the origin of these residuals. 
A disc reflection model with the Galactic cold absorption 
column fixed at zero (rather than 1.55 x 10 20 cm -2 , from 
Table can account for the shape of these residuals, with 
a xl °f 1-066 for 778 degrees of freedom. A simple model fit 
with no Galactic absorption is still inadequate, with a xt 
of 1.264 for 782 degrees of freedom. These fits are unlikely 
to be physical, the Galactic column is rarely lower than 
~ 5 x 10 19 cm -2 . The residuals are instead likely to be due 
to some curvature of the soft excess not addressed by our 
models. 
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Table 3. Rclativistically-blurrcd photoionizcd disc fits to the sources. Em. Ind. is the index of the emissivity of the accretion disc 
(assumed to be a power law), Ri n is the disc's inner radius (in gravitation radii), i is the inclination of the disc to the line of sight (in 
degrees), Fe is the iron abundance of the disc (relative to solar), § (in erg cm s —1 ) is the ionization parameter of the gas (£ = 4nF/n[i, 
where F is the illuminating energy flux and nu is the hydrogen number density in the illuminated layer). Flux Frac. is the flux fraction, 
see Section HI 
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dominated and includes no power law component. c Adding non-redshifted (local to our Galaxy) edges improves the fit, see Section |3 
The model based on a non-rotating black hole has a similar goodness of fit, see Table e The fit includes warm absorption edges or a 

narrow Gaussian line at 6.4 keV, with parameters listed in Table HI 

Table 4. Redshifted absorption edges and narrow 6.4 keV Gaussian emission lines included in the relativistically-blurred disc fits in 
Table |3 Energy is the energy of the absorption edge, r is the optical depth and Width is the equivalent width of the Gaussian line. 
Energies given without associated errors are fixed at 0.74 or 0.87 keV, the energy of the O VII and O VIII K absorption edges. Blank 
spaces indicate no feature is present. 
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Table 5. Photoionized disc fits to the sources where the relativistic blurring is calculated assuming a non-rotating black hole. Column 
labels are as Table El 
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Figure 4. The residuals to the simple model (left) and relativistically-blurred photoionized disc reflection model (right) fits to PG 
0050+124, PG 0844+349 and PG 1115+407. The disc reflection model fit to PG 1115+407 does not include a power law component. 
The fit parameters are given in Tables 121 151 and HI The disc fits to these sources are an improvement on the simple fits, with xi (degrees 
of freedom) going from 1.270 (928) to 1.117 (924), 1.041 (631) to 0.987 (627) and 0.963 (416) to 0.913 (413) for PG 0050+124, PG 
0844+349 and PG 1115+407, respectively. 



5.3 PG 1404+226 

PG 1404+226 has been previously fit with the disc reflection 
model in Crummy et al. (2005). PG 1404+226 is strongly 
variable, and when the simple model is used it appears to 
have an absorption feature in the low state which disappears 



in the high state (see also Dasgupta et al. 2005). Crummy 
et al. show that the disc reflection model can reproduce 
both the high and low spectral states without the need for 
variable absorption. 
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Figure 5. The relativistically-blurred ionized disc reflection model plotted in vF v , i.e. showing the flux emitted as a function of energy. 
The black crosses are the data points, the blue line is the fit to the data (from Table [3] with the absorption edges for PG 1309+355 
(0.74 keV) and ARK 564 (0.68 keV) from Table H"! . the green line is the power law component of the fit, the red line is the reflection 
component of the fit, and the insert shows the Laor profile (also plotted as vF v ) that the two components are blurred with to produce 
the model. The figure illustrates the effect of ionization parameter (£) on the model, the model shown for PG 1309+355 has f = 3, PG 
1501+106 has £ = 510, PG 1116+215 has £ = 1270 and ARK 564 has £ = 3120. The lines are clearly more distinct in the reflection 
component at lower ionization parameters, and the Compton reflection continuum increases as ionization parameter increases. An iron 
line between 6.4 and 6.7 keV is clear in the reflection component for all values of the ionization parameter, but can only be seen in the 
source spectra as a slight enhancement to the continuum as it is blurred by the wide Laor lines (see Gallo et al. 2004b for an investigation 
of this feature in MRK 0586). 



5.4 NGC 4051 

The observation of NGC 4051 we analyse is of very high 
quality, and shows some features not observable in any 
of the other sources. It is clear from Fig. |7] that the disc 
reflection model follows the shape of the continuum very 
well, but there are several small features which contribute 
to the fairly high xt of 1.240 for 1040 degrees of freedom. 
Some of these features may be calibration problems, such 
as around the XMM-Newton instrumental Au M edge at 
2.3 keV, and some may be due to small amounts of absorp- 
tion or emission. However it is also possible that some of the 
approximations made in the disc reflection model are invalid 
when analysing observations of this detail, see Section 
A more detailed analysis of NGC 4051 and its variabil- 
ity using the disc reflection model is in Ponti et al. (in prep.). 



5.5 E 1346+266 

E 1346+266 is unusual in our sample as it has a high 
redshift, z = 0.915. This enables us to investigate the 
effectiveness of the disc reflection model at source energies 



up to 23 keV. The simple model, the disc reflection model 
and the version of the disc reflection model with a non- 
rotating black hole are all roughly equal in quality of fit to 
the source, so whilst the disc reflection model has proved 
itself to work at least as well as the simple model in this 
energy range it is difficult to say anything definite about 
the source. The disc reflection model essentially works by 
fitting the soft excess and any excess iron emission. The soft 
excess is less visible in high redshift AGN due to being both 
shifted out of the instrument band and strongly absorbed. 
The disc reflection model will therefore be most useful in 
high redshift sources which have excesses in the continuum 
around 6.4 keV due to iron, although it is able to fit sources 
where this is not the case. 



5.6 PG 1211+143 

PG 1211+143 has a feature in the 4-8 keV region which 
is not predicted by either model, see Fig. |H| The shape 
of the residuals to the fits suggest two possibilities, an 
absorption edge at ~7 keV or a broad Laor iron line. We 
investigated the latter possibility by adding a Laor line 
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Figure 7. The simple model (left) and relativistically-blurred photoionized disc reflection model (right) fits and residuals to NGC 405f , 
including Gaussian lines and absorption edges as stated in Tables [2] [3] and [3} . The disc reflection model is clearly a superior fit, and 
follows the shape of the continuum very well, despite the high xt- The spectra have been re-binned to show the shape of the residuals 
more clearly. 
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Figure 6. A histogram showing the measured flux fractions of 
the sources. All of the sources have a reflection component, and 
many sources are dominated by reflection. The deficit of sources 
in the range 0.8 — 0.9 may be partly due to low data quality, see 
Section [I] 



at 6.4 keV (in the source frame) to our reflection model, 
with the inclination, emissivity index and outer radius 
fixed to the same values as our Laor blurring, allowing 
only the inner radius and normalisation to vary. We found 
this improved the fit very significantly, to a xt of 1.053 for 
932 degrees of freedom. However, this is unphysical as any 
iron line would also have other reflected lines associated 
with it. We therefore fit PG 1211+143 with two reflection 
components, both illuminated by the same power law and 
at the same inclination, with the same iron abundance, 
but allowing all other parameters to vary. We found that 
this improves the fit further, to a xi °f 1-030 for 930 
degrees of freedom, with illuminating power law index 
2-07 -o' oi, non-redshifted edges (i.e. potentially local to our 
Galaxy) at 0.74 and 0.87 keV with depths of 0.28±° ° l and 



0.241q °4 respectively, disc inclination 47.6^04 degrees, iron 
abundance 6.4j^'2 times solar, emissivity indices 9.21q3 
and 4.8lJ'g, inner radii 1.311q'q7 and 3.01q'7 gravitational 
radii, and ionization parameters 6IOI4Q and 29lj erg 
cm s _1 . Where we give two parameters the second refers 
to the outer disc component, with an inner radius of 3 
gravitational radii. 

The fit parameters are not consistent with the fit in Table 
[3] but given that the fit has improved so much (Axl of 
235 for 7 degrees of freedom) that is entirely expected; in 
particular the iron abundance has reduced from 10 to 5 
since a single ionized disc is no longer trying to fit both the 
large iron line and the soft excess. This model is somewhat 
similar to the work of Sobolewska & Done (2004), who 
achieve a similar goodness of fit to PG 1211+143 with 
three unblurred reflection components. Taking account of 
relativity by adding relativistic blurring removes the need 
for one of their reflectors, we instead have a power law. Our 
model is also close to that of Fabian et al. (2002) who use 
three blurred reflectors to model 1H 0707-495. 
We also investigated the alternative possibility already re- 
ported by Pounds et al. (2003) which is an absorption edge 
at ~7 keV. We added a redshifted edge at 7 keV with energy 
and depth as free parameters to our reflection model, and 
found that it too improved the fit significantly. This model 
has a °f 1-122 for 934 degrees of freedom, with emissivity 
index 6.48l l ' o l 5, inner radius 1.24io'oo gravitational radii, 
inclination 35.31^4 degrees, iron abundance 4.73±^ times 
solar, illuminating spectral index T881q'qii ionization 
parameter 3OOI3 erg cm s _1 , edge energy 7.44l l ' l j'o keV 
and edge depth 0.84lg'Qg. A plausible way to create the 
absorption edge is a highly ionized outflow (Pounds et al. 
2003). We find that adding an edge at ~7.3 keV to the 
disc reflection model is an improvement, but that adding 
a second reflection component at 6.4 keV gives an even 
better fit. Physically, since this second component has a 
larger inner radius, it may represent a a hot spot or ring of 
emission excited by e.g. flares from the magnetic fields. 
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5.7 PHL 1092 

PHL 1092, when fit with the disc reflection model, appears 
to have a non-redshifted edge at 0.87 keV. The fit param- 
eters when this edge is included are a xt °f 1-208 for 179 
degrees of freedom (compared to 1.237 for 180 degrees of 
freedom without the edge), edge depth 0.31q2, emissivity 
index 9.£T_2'2> inner radius 1.24+q'q 8 gravitational radii, in- 
clination 69.4^g 2 degrees, iron abundance l.OOj^os 

times 

solar, illuminating spectral index 2.53_Zoloi an d ionization 
parameter 1460t 3 Q er S cm s_1 - 

PHL 1092 has been fit by Gallo et al. (2004a), who use 
the radiative efficiency formula of Fabian (1979) to conclude 
that the variability of the source is inconsistent with a non- 
rotating black hole, whilst being consistent with a maximally 
rotating black hole. They suggest that the source is reflec- 
tion dominated, and that light bending may be in operation. 
This is consistent with our best fit for PHL 1092, which is 
pure reflection without a power law component. 



5.8 MRK 0586 

MRK 0586, when fit with the disc reflection model, appears 
to have a non-redshifted edge at 0.87 keV. The fit param- 
eters when this edge is included are a xt °f 0.977 for 522 
degrees of freedom (compared to 1.007 for 523 degrees of 
freedom without the edge), edge depth 0.151q'o5i emissivity 
index 4.8jl5'i, inner radius 1.24+q'qq gravitational radii, 
inclination 51.3lg Q degrees, iron abundance 0.75_to'io times 
solar, illuminating spectral index 2.33 jig' 02 an d ionization 
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1 2 5 
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parameter 1710l2oo er S cm s_i - 
MRK 0586 has been previously investigated by Gallo et 
al. (2004b) who find a feature consistent with a highly 
broadened iron line emitted from the inner regions of an 
accretion disc around a Kerr black hole. Our disc fit, where 
the fit is dominated 3 by the soft excess rather than the weak 
iron line, supports their hypothesis; we find an inclination 
(55l^ degrees) consistent with theirs (37_^ 4 ) although we 
have fit the emissivity index and find 4.2IJ 4 where they fix 
it at 3. 
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Figure 8. From top to bottom the simple model, disc reflection 
model and double relection disc reflection model fits and residuals 
to PG 1211+143, and a plot of the double reflection disc reflection 
model (the solid line is the total model, the dotted line is the 
power law and the dashed and dash-dotted lines are the reflection 
components). PG 1211+143 has some structure in the 4-8 keV 
region, visible in the top two spectra, that is best explained by 
adding a second reflection component with an inner radius of 3 
gravitational radii, as shown in the bottom panels. See Section IS1 
for a discussion of this feature. The spectra have been re-binned 



6 DISCUSSION 

The relativistically-blurred photoionized disc reflection 
model is clearly a better fit in terms of reduced chi-squared. 
25 of the 34 sources show a significant (Axl > 2.7 per 
degree of freedom; note that this does not correspond to 
a 90 per cent probability as the x 2 distribution is not 
calibrated between models) improvement when moving 
from the simple model to the disc reflection model, with 6 
sources showing a significant worsening and the remaining 
3 sources inconclusive. The disc reflection model is a good 
fit to the sources (xt < 1-1) in 22 of the sources, an 
adequate fit (xt < 1.3) in 10 cases, and a poor fit (xt > 
1.3) to TON S180, compared to 19 good fits, 13 adequate 
fits and 2 poor fits with the simple model. The sources 



3 The majority of the photon flux is in the soft energy band, 
and fitting just the soft band recovers parameters consistent with 
fitting over the entire available energy range. 
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which are fit worse by the disc reflection model are PG 
1211+143, PG 1309+355, IRAS 13349+2438, PHL 1092, 
RE J1034+396 and TON S180. PG 1211+143 can be 
well fit by the disc reflection model if a second reflection 
component is added (Section the other sources show no 
obvious commonalities, being both high and low quality 
observations with fit parameters spread out over a similar 
range to the entire sample. Several of these sources show 
small excesses around 6 keV with a drop at 7 - 8 keV, 
which may be due to extra iron emission like PG 1211+143, 
however PG 1309+355 does not, and this feature is not 
uncommon throughout the entire sample. 
The disc reflection model accounts naturally for the shape 
of the continuum in all the sources. The model also natu- 
rally accounts for features which would be interpreted as 
absorption edges in the simple model - 17 sources appear 
to have absorption edges when the simple model is used, 
reducing to 7 sources using the disc reflection model. We 
therefore conclude that the relativistically-blurred pho- 
toionized disc reflection model reproduces the spectra of our 
sources better than the simple model, and should be a pri- 
mary tool in the investigation of X-ray observations of AGN. 



6.1 Warm absorption 

Only a few sources may show warm absorption edges local 
to our Galaxy. When using the simple model we find only 
one source where a non-redshifted (i.e. local) edge improves 
the fit, and we eliminate that as a possibility based on a 
variability analysis in a previous paper (see PG 1404+226 
in Section |SJ. When using the disc reflection model we 
find three sources that may exhibit non-redshifted edges 
(PG 1211+143, PHL 1092 and MRK 0586), none of which 
have warm absorption edges local to the sources. Whilst 
our analysis is not optimised to detecting these features, 
we can conclude that strong non-redshifted edges are rare. 
A further point to note is that the disc reflection model 
finds more of them than the simple model, which is the 
opposite of the situation with edges local to the sources. 
Warm absorption edges local to the source appear to be 
very common when the simple model is used, but many of 
the features that seem to be edges can be explained by the 
shape of the disc reflection model. Many features previously 
identified as absorption edges may instead be due to blurred 
reflection from an ionized disc; high resolution observations 
and variability studies offer methods to distinguish between 
the two cases. 



6.2 Inclinations 

The model parameters we measure give us information 
about our sample. The ranges of inclinations are plotted 
in Fig. [5] It should be noted that these inclinations are 
rather model dependent, and that the formal errors from 
xspec fitting are likely to be underestimates. The inclina- 
tion of the accretion disc to the line of sight spans the a 
range from 22° - 90°, though three sources have a incli- 
nations consistent with 0° within the errors (see Table 0. 
A random distribution in inclination would be proportional 



to sin(i), which is not consistent with our measurements. 
A Kolmogorov-Smirnov (K-S) test measures the probability 
that a set of data is consistent with a given distribution. A 
K-S test against a random distribution to our data, neglect- 
ing the measured uncertainties in inclination, gives a prob- 
ability of 0.06. Under the unified model of active galactic 
nuclei (Antonucci 1985) they would be confined to a range 
in inclination, P(i) oc sin(i); i < a where 90 — a is the angu- 
lar extent of the torus above the plane of the disc in degrees. 
The unified model can explain the apparent deficit in discs 
with an inclination above 70°, though it does not explain 
the sources which we do measure with high inclinations, nor 
the apparent deficit below 30° . We performed a series of K- 
S tests (in OCTAVE 2.1.60) using our simple unified model 
distribution with cut-offs a ranging from 0° to 50°, i.e. over 
a range in maximum inclinations of 40° - 90°, discarding 
any measured inclinations above the cut-off. The results are 
shown in Fig. 1101 We found a maximum K-S test probabil- 
ity of ~ 0.34 using an 81° cut-off, and probabilities of ~ 0.1 
with cut-offs over the range 50° - 85°. The data are incom- 
patible with being randomly distributed over the entire 0° - 
90° range in inclination, but are somewhat consistent with 
being random over 0° - 80°. The deficit below 30° may be 
partly explained by a selection effect; under the relativistic 
disc model sources emit most of their light along the disc 
and very little directly upwards, due to Doppler beaming. 
It is also possible that it is an artefact of our small sample 
size. 

The emission from the relativistically-blurred photoionized 
disc reflection model is dominated by the innermost part of 
the accretion disc, if the inner accretion disc is not always 
coplanar with the torus (Bardeen & Petterson 1975) then 
our high inclination sources can be explained; however, esti- 
mates of the time taken for a Bardeen-Petterson inner disc 
to align with the outer disc are 10 6 years (Natarajan & 
Armitage 1999) so Bardeen-Petterson discs should be rare. 
It is also possible that some of these type 1 AGN have a 
small or no torus. 



6.3 Black hole rotation 

For every source in our sample we found that using rela- 
tivistic blurring based on a non-rotating black hole was a 
worse fit than a maximally rotating hole (Section [3J, with 
only two sources where the fit is not significantly worse. 
This implies that rapidly rotating black holes dominate our 
sample. More information can be extracted by considering 
the measured inner radii of our sources. An accretion disc 
around a black hole is not expected to extend further in 
than the last stable orbit. The last stable orbit around a 
non-rotating (Schwarzschild) black hole is at 6 gravitational 
radii. All 34 sources have inner accretion disc radii consis- 
tent with being under 6 gravitational radii, and only two 
sources (PG 1244+026 and E 1346+266, the same sources 
that are equally consistent with a non-rotating fit) have an 
inner radius which may be over that limit. The last stable 
orbit around a maximally rotating (Kerr) black hole is at 
1.235 gravitational radii (this requires a disc with angular 
momentum in the same direction as the black hole angular 
momentum). Almost all of the sources (29 of the 34) have 
inner radii consistent with being below 1.3 gravitational 
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Figure 9. A histogram showing the measured inclinations of 
the sources (solid line) with the random distribution over-plotted 
(dotted line). There is a deficit of sources at high inclinations 
compared to the random distribution, which supports the AGN 
unified model. 



radii. This shows that the majority of black holes in (type 
1) AGN are rapidly rotating, in agreement with theoretical 
predictions based on the fact that black holes gain most 
of their mass from accretion discs which tend to increase 
the angular momentum. Volonteri et al. (2005) predict 
a distribution of spins where 70 per cent of AGN black 
holes are maximally rotating, which agrees well with our 
measurements. 

Black hole rotation has also been measured for more distant 
quasars by Streblyanska et al. (2005), who find an iron 
line described by a Laor profile with an inner radius of 
~3 gravitational radii in a combined spectrum of distant 
quasars from the Lockman Hole. This implies that the 
average distant AGN rotates rapidly. 

It should be noted that the black hole rotation is not a free 
parameter in our analysis; we perform the fits assuming 
a maximally rotating black hole therefore we have not 
measured the rotation, merely found that many of the 
fits are statistically consistent (low xi) with maximal 
rotation. An alternative possibility is that radiation from 
inside the last stable orbit (the plunging region) may be 
important. See Krolik & Hawley (2002) for a discussion of 
the inner edges of accretion discs and Section 3.4 in Fabian 
& Miniutti (2005) for comments on why the inner radius is 
likely to be a good estimator of black hole spin. 
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Figure 10. The Kolmogorov-Smirnov probability (Pk—s) that 
the measured inclinations are consistent with a random distri- 
bution over the range 0° to a maximum inclination (x-axis). The 
data are inconsistent with being random over the entire 90° range, 
though they are reasonably consistent with being random over 0° 
— 80°, giving support to an AGN unified model with a with a 
torus that obscures discs of inclinations above 80° . 



6.4 Emissivity indices 

The emissivity index of our sources tends to be very high, 
with 22 sources having an index over 6, 11 sources with an 
index between 3 and 6, and only one, E 1346+266, with an 
index below 3. The measurements for E 1346+266 are not 
strongly constrained as it is at a fairly high redshift (see 
Section and the data are consistent with a much higher 
index. For a typical source with an emissivity index of 7, 
inclination of 50°, and inner radius of 1.235 gravitational 
radii an observer will detect the radial flux profile shown in 
Fig. 1111 Energy release is very concentrated towards the cen- 
tre; 90 per cent of the flux is emitted within 2.0 gravitational 
radii. The energy observed is slightly less concentrated, with 
90 per cent of the flux originating within 3.2 gravitational 
radii. This is different to the result of naively integrating the 
emissivity profile due to relativistic effects, however it is still 
clear that we are probing very close to the central black hole. 
This suggests that a study to compare measured variability 
timescales with emission region sizes from the disc reflection 
model could prove fruitful. 



6.5 Abundances 

The iron abundance measured from the sources is consis- 
tent with being between solar and 1/3 solar for 28 of the 
34 sources. Six sources have iron abundances which are 
definitely above solar. These fairly low abundances partially 
explain the lack of obvious broad iron lines. Care should be 
taken in drawing strong conclusions about metallicities from 
this as the model only allows for a varying iron abundance, 
the abundances of all other elements are fixed at solar. 
Furthermore, fits with the iron abundance fixed at a high 
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Figure 11. A cumulative plot of the flux in a Laor line with 
parameters given in Section |§] emitted as a function of radius. 
This was generated by varying the inner radius parameter and 
normalising to the flux emitted when the inner radius is 1.235 
gravitational radii. The small step-like features in the plot are 
due to numerical inaccuracies in the model. 
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Figure 12. A plot comparing spectral indices measured using 
both models. r s ; mp ; e is the spectral index from the simple model, 
Tdisc is the spectral index from the relativistically-blurred pho- 
toionized disc reflection model. The dotted line shows a 1:1 ratio. 
Most sources have a similar spectral index under both models, 
but the simple model has several outliers with a lower spectral 
index than expected. 



value are still acceptable in some cases. However, these 
results encouragingly show that this technique is capable of 
measuring metallicities of accretion discs, and with higher 
quality data and a model which allows the measurement 
of abundances for more elements the question of the 
composition of AGN accretion discs could be answered by 
observing their X-ray emission. 



6.6 Spectral indices 

The spectral index we measure from the sources using both 
models is shown in Fig. 1121 Most of the sources have simi- 
lar indices under both models, but there are several sources 
which have much lower spectral indices under the simple 
model (PG 1202+281, PG 1307+085, PG 1404+226, NGC 
4051, and PHL 1092). PHL 1092 has no power law compo- 
nent in its spectra when fit with the disc reflection model, 
and a very strong soft excess, so the disc reflection model 
has a steeper illuminating continuum than the simple model. 
The other four sources have reflection components which in- 
crease in strength above ~8 keV so the simple model, which 
effectively fits the entire high energy region with just a power 
law, will be flatter. See Fig. 1131 These results show that cau- 
tion should be exercised when measuring the spectral index 
of a source if the intention is to relate it to a physical model; 
if disc reflection is present the spectral index will be due to 
the sum of more than one processes. 
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Figure 13. The relativistically-blurred photoionized disc reflec- 
tion model used to fit PG 1202+281, shown over an extended 
energy range. The dashed line is the reflection component, the 
dash-dotted line is the illuminating power law, and the solid line 
is the total model. The curvature at low energies is due to Galac- 
tic cold absorption. See Table l3l for fit parameters. Note that the 
slope of the overall spectrum is significantly different from the 
slope of the power law component due to the high energy upturn 
in the reflection component (see Spectral indices in Section ISV 



6.7 The relativistically-blurred photoionized disc 
reflection model 

The relativistically-blurred photoionized disc reflection 
model is an improvement in both goodness of fit and 
physical motivation over the simple model. However, in its 
present form there are still residuals to the fit in several 
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of the sources. These residuals tend to occur in the 6 - 
8 keV region, and show a variety of features: some may be 
absorption edges or notches, some emission lines (narrow 
and broad). The disc reflection model is very smooth due to 
the relativistic blurring, it cannot produce a notch, though 
the model can explain features that look like absorption 
edges. The sharper features are therefore probably due 
to iron (and possibly other elements like nickel) emission 
and absorption, with the broader features due to either 
absorption from gas in the line of sight or iron emission 
from further out in the disc - perhaps from regions excited 
by magnetic reconnection. It is possible that our current 
model of a geometrically flat accretion disc with a power 
law emissivity index does not accurately represent all the 
AGN studied. 

The photoionized reflection spectra are generated by 
a model of illumination onto a slab of gas of constant 
density. This is unlikely to be the situation in a real 
AGN accretion disc, but is a necessary approximation to 
allow the calculation of a large grid of models such that 
fits can be performed with xspec. Previous studies have 
compared constant density reflection spectra to those which 
use a more physical density distribution (e.g. hydrostatic 
equilibrium, which still requires assumptions about the 
structure of the disc and other factors like the external 
coronal pressure). Ballantyne, Ross & Fabian (2001) found 
that their hydrostatic equilibrium model could be approx- 
imated well by a "diluted" constant density model, where 
diluted means that the reflection fraction of the constant 
density model is smaller than the reflection fraction of 
the hydrostatic model it is fitting. This happens because 
an outer layer of the disc becomes highly ionized (a low 
density means a high ionization parameter) to the point 
where it produces no line features, and the lines form in a 
partially hidden layer below this. This may mean that the 
flux fractions we measure in Section ^. fl are too low, some of 
the emission we interpret as due to the illuminating power 
law might be reflection from a highly ionized surface layer. 
An accretion disc is an uncontroversial component in a 
model of an AGN, but our model involves an accretion 
disc illuminated by a power law. This raises the question 
of the origin of the illuminating radiation. A corona of hot 
electrons above the accretion disc is often hypothesised, 
and can produce a power law by Comptonisation of UV 
disc photons (or bremsstrahlung or synchrotron emission). 
Alternatively a small hot-spot, either the base of a jet, 
shocks in a failed jet (Ghisellini, Haardt & Matt 2004) or 
a region excited by a magnetic flare above the disc, can 
produce the necessary radiation. This model seems more 
likely in the context of our results - we have several sources 
where the illuminating continuum is not visible to us (Table 
[3] Figure [fjj and this can be explained by the light bending 
effects of the central black hole - photons emitted from 
sources above the disc tend to be bent towards it and very 
few escape to be observed, increasing the measured flux 
fraction by supressing the power law. Reflected radiation 
escapes the black hole more easily, it is Dopper beamed 
along the plane of the disc, away from the black hole. This 
is discussed in Miniutti et al. (2003) and Miniutti & Fabian 
(2004). However, it should be noted that our results show 
no correlation between flux fraction and inclination, so 
we consider that Doppler beaming along the disc is not 



the dominant factor. This would also be the case if the 
illuminating hot spot rotates with the disc, sharing the same 
velocity and beaming; unlikely for jet based illumination 
but expected for magnetic flares above a disc or a corona. 
The consistently small inner radius is difficult to explain 
with magnetic flares. Alternatively, with a small power law 
emitting region it is possible that the power law radiation 
is obscured, although since I in 3 of the sources have no 
power law component and all sources have disc components 
the geometry of the absorber is difficult to imagine. Our 
favoured scenario is a small power law illuminating region 
above the central black hole, where light bending effects 
mean that the illuminating power law radiation can be 
deflected onto the disc and not be detected.. 
One further problem that the relativistically-blurred pho- 
toionized disc reflection model can solve is the general 
absence of strong iron lines in AGN spectra. An accretion 
disc illuminated by X-rays should produce an iron line. The 
disc reflection model shows that the iron line is present, but 
relativistically blurred to such an extent that it is hard to 
detect, see Figure |S] 



7 CONCLUSIONS 

We investigated a large sample of type 1 AGN X-ray spectra, 
and fit them with the relativistically-blurred photoionized 
disc reflection model of Ross & Fabian (2005). This work 
has shown that: 

• The relativistically-blurred photoionized disc reflection 
model fits the sample better than the conventional simple 
model. Since the disc reflection model is physically moti- 
vated, unlike the simple model, it can also give us direct 
information about the inner accretion disc. 

• The disc reflection model reproduces the continuum 
shape of all the AGN in the sample, in particular the soft 
excess is naturally explained by the model. The soft excess 
is composed of many broad lines mixed together. 

• The disc reflection model reproduces many features in 
the spectrum that could otherwise be interpreted as warm 
absorption edges. 

• The majority of black holes in the sample are strongly 
rotating, alternatively it is possible that radiation from the 
plunging region is important. 

• The inclinations measured from the sample span the 
range of 20° - 90°, so neither edge of the distribution is 
consistent with the simplest unified model. However, there 
is a deficit at inclinations > 70° so there is some evidence 
for torus obscuration. 

• The disc reflection model provides evidence that the 
iron abundances of black hole accretion discs tend to be 
solar or mildly sub-solar. 

We have shown that the relativistically-blurred photoion- 
ized disc reflection model is an important tool in the study 
of AGN, and that taking account of the intrinsically rela- 
tivistic nature of these sources not only provides information 
about the state of the system, but explains the apparent ab- 
sence of iron line reflection from the disc and the shape and 
behaviour of the soft excess. 
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